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Abstract We discuss suitable classes of diffusion processes, for which functionals 
relevant to finance can be computed via Monte Carlo methods. In particular, we 
construct exact simulation schemes for processes from this class. However, should 
the finance problem under consideration require e.g. continuous monitoring of the 
processes, the simulation algorithm can easily be embedded in a multilevel Monte 
Carlo scheme. We choose to introduce the finance problems under the benchmark 
approach, and find that this approach allows us to exploit conveniently the analytical 
tractability of these diffusion processes. 



1 Introduction 

In mathematical finance, the pricing of financial derivatives can under suitable con- 
ditions be shown to amount to the computation of an expected value, see e.g. [53], 
[56]. Depending on the financial derivative and the model under consideration, it 
might not be possible to compute the expected value explicitly, however, numerical 
methods have to be invoked. A candidate for the computation of such expectations 
is the Monte Carlo method, see e.g. [11], [30], and [44]. Applying the Monte Carlo 
method typically entails the sampling of the distribution of the relevant financial 
state variables, e.g. an equity index, a short rate, or a commodity price. It is then, 
of course, desirable to have at one's disposal a recipe for drawing samples from the 
relevant distributions. In case these distributions are known, one refers to exact sim- 
ulation schemes, see e.g. [55], but also [7], [8], [9], and [16], for further references 
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on exact simulation schemes. If exact simulation schemes are not applicable, dis- 
crete time approximations, as analyzed in [44] and [55] become relevant. In recent 
years, it has been shown under certain assumptions that using the multilevel Monte 
Carlo method, see [29] and also [38], [39], the standard Monte Carlo convergence 
rate, achieved by exact simulation schemes, can be recovered. 

For modeling financial quantities of interest, it is important to know a priori if 
exact simulation schemes exist, so that financial derivatives can be priced, even if 
expected values cannot be computed explicitly. In this paper, we discuss classes of 
stochastic processes for which this is the case. For one-dimensional diffusions, Lie 
symmetry analysis, see [10], and [54] turns out to be a useful tool. Besides allowing 
one to discover transition densities, see [21], it also allows us to compute Laplace 
transforms of important multidimensional functionals, see e.g. [20]. In particular, 
we find that squared Bessel processes fall into the class of diffusions that can be 
handled well via Lie symmetry methods. 

The Wishart process, [13], is the multidimensional extension of the squared 
Bessel process. It turns out, see [33] and [34], that Wishart processes are affine pro- 
cesses, i.e. their characteristic function is exponentially affine in the state variables. 
We point out that in [33], and [34] the concept of an affine process was generalized 
from real-valued processes to matrix-valued processes, where the latter category 
covers Wishart processes. Furthermore, the characteristic function can be computed 
explicitly, see [33], and [34]. Finally, we remark that in [1] an exact simulation 
scheme for Wishart processes was presented. 

Modeling financial quantities, one aims for models which provide an accurate 
reflection of reality, whilst at the same time retaining analytical tractability. The 
benchmark approach, see [56], offers a unified framework to derivative pricing, risk 
management, and portfolio optimization. It allows us to use a much wider range of 
empirically supported models than under the classical no-arbitrage approach. At the 
heart of the benchmark approach sits the growth optimal portfolio (GOP). It is the 
portfolio which maximizes expected log-utility from terminal wealth. In particular, 
the benchmark approach uses the GOP as numeraire and the real world probability 
for taking expectations. We find that the class of processes for which exact simula- 
tion is possible is easily accommodated under the benchmark approach, which we 
illustrate using examples. 

The remaining structure of the paper is as follows: In Section 2 we introduce the 
benchmark approach using a particular model for illustration, the minimal market 
model (MMM), see [56]. Section 3 introduces Lie symmetry methods and discusses 
how they can be used in the context of the benchmark approach. Section 4 presents 
Wishart processes and shows how they can be used to extend the MMM. Section 6 
concludes the paper. 
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2 Benchmark Approach 

The GOP plays a pivotal role as benchmark and numeraire under the benchmark 
approach. It also enjoys a prominent position in the finance literature, see [43], but 
also [12], [45], [42], [49], [50], and [58]. The benchmark approach uses the GOP 
as the numeraire. Since the GOP is the numeraire portfolio, see [49], contingent 
claims are priced under the real world probability measure. This avoids the restric- 
tive assumption on the existence of an equivalent risk-neutral probability measure. 
We remark, it is argued in [56] that the existence of such a measure may not be 
a realistic assumption. Finally, we emphasize that the benchmark approach can be 
seen as a generalization of risk-neutral pricing, as well as other pricing approaches, 
such as actuarial pricing, see [56]. 

To fix ideas in a simple manner, we model a well-diversified index, which we 
interpret as the GOP, using the stylized version of the MMM, see [56]. Though 
parsimonious, this model is able to capture important empirical characteristics of 
well-diversified indices. It has subsequently been extended in several ways, see e.g. 
[56], and also [4]. To be precise, consider a filtered probability space (Q,£?,£/_,P), 
where the filtration s[_ = (M)te[o,°°) i s assumed to satisfy the usual conditions, 
which carries for simplicity one source of uncertainty, a standard Brownian mo- 
tion W = {W t , t G [0,°°)}. The deterministic savings account is modeled using the 
differential equation 

dS? = rS°dt, 

for t G [0,oo) with Sq = 1, where r denotes the constant short rate. Next, we intro- 
duce the model for the well diversified index, the GOP S t \ which is given by the 
expression 

S f* =S 0s?*=S°Y t a t S *. (1) 

a?* 

Here Y t = is a square-root process of dimension four, satisfying the stochastic 
differential equation (SDE) 

dY t = (l-rjY t )dt + y/Y t dW t , (2) 

for t G [0,°o) with initial value Yq > and net growth rate 77 > 0. The deterministic 
function of time oc, * is given by the exponential function 

a, 5 * = aoexp{T7f} , 

with scaling parameter Oq > 0. Furthermore, it can be shown by the Ito formula that 
a, * is the drift at time t of the discounted GOP 
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so that the parameters of the model are Sq' , Oq, r\, and r. We note that one obtains 
for the GOP the SDE 

dsf* = S?* Ur+^)dt+ \[^dW^j , (3) 

which illustrates the well-observed leverage effect, since as the index S, * decreases, 
its volatility -j= = <J J^- increases and vice versa. 
It is useful to define the transformed time <p(t) as 

<P(t) = <p(0) + \J o af'ds. 

Setting 

we obtain the SDE 

dX 9(f) = 4d(p{t) + 2^X^~ ) dW (p{t) , (4) 

where 

dw v(t) = y -j-dWr 

for t £ [0,o°). This shows that X = {Xy, (p £ [<p(0),°°)} is a time transformed 
squared Bessel process of dimension four and W = {Wy, (p £ [<p(0),°°)} is a Wiener 
process in the transformed (p-time (p(t) £ [<p(0),°°), see [57]. The merit of the dy- 
namics given by (4) is that transition densities of squared Bessel processes are well 
studied; in fact we derive them in Section 3 using Lie symmetry methods. 

We remark that the MMM does not admit a risk-neutral probability measure be- 

s s * 

cause the Radon-Nikodym derivative A t = of the putative risk-neutral measure, 

$t 

which is the inverse of a time transformed squared Bessel process of dimension 
four, is a strict local martingale and not a martingale, see [57]. On the other hand, 
S s *, is the numeraire portfolio, and thus, when used as numeraire to denominate 
any nonnegative portfolio, yields a supermartingale under the real-world probabil- 
ity measure P. This implies that the financial market under consideration is free 
of those arbitrage opportunities that are economically meaningful in the sense that 
they would allow to create strictly positive wealth out of zero initial wealth via a 
nonnegative portfolio, that is, under limited liability, see [48] and [56]. This also 
means that we can price contingent claims under P employing S s * as the numeraire. 
This pricing concept is referred to as real-world pricing, which we now recall, see 
[56]: for a nonnegative contingent claim with payoff H at maturity T, where H is 

sifr -measurable, and E ^^r^j < °° > we define the value process at time t £ [0, T] by 
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V t := sf'E 




) 



(5) 



Note that since Vt = H, the benchmarked price process is an (^,f)-martingale. 



Formula (5) represents the real-world pricing formula, which provides the minimal 
possible price and will be used in this paper to price derivatives. If the expectation in 
equation (5) cannot be computed explicitly, one can resort to Monte Carlo methods. 
In that case, it is particularly convenient, if the relevant financial quantities, such as 
S T * can be simulated exactly. In the next section, we derive the transition density of 
S * via Lie symmetry methods, which then allows us to simulate S T * exactly. Note, 
in Section 4, we generalize the MMM to a multidimensional setting and present a 
suitable exact simulation algorithm. 



3 Lie Symmetry Methods 

The aim of this section is to present Lie symmetry methods as an effective tool 
for designing tractable models in mathematical finance. Tractable models are, in 
particular, useful for the evaluation of derivatives and risk measures in mathematical 
finance. We point out that in the literature, Lie symmetry methods have been used to 
solve mathematical finance problems explicitly, see e.g. [19], and [40]. Within the 
current paper we want to demonstrate that they can also be used to design efficient 
Monte Carlo algorithms for complex multidimensional functionals. 

The advantage of the use of Lie symmetry methods is that it is straightforward 
to check whether the method is applicable or not. If the method is applicable, then 
the relevant solution or its Laplace transform has usually already been obtained in 
the literature or can be systematically derived. We will demonstrate this in finance 
applications using the benchmark approach for pricing. 

We now follow [20], and recall that if the solution of the Cauchy problem 



is unique, then by using the Feynman-Kac formula it is given by the expectation 



where Xq = x, and the stochastic process X = {X t , t > 0} satisfies the SDE 



u t = bx^Uxx + f(x)u x — g{x)u , x > , t > , 
u(x,0) = (p(x),x£ Q = [0,°°), 



(6) 
(7) 
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We now briefly describe the intuition behind the application of Lie Symmetry meth- 
ods to problems from mathematical finance, in particular, the integral transform 
method developed in [47], and the types of results this approach can produce. Lie's 
method allows us to find vector fields 

v = %(x,y,u)dx + T(x,t,u)dt + (j)(x,t ,u)d u , 

which generate one parameter Lie groups that preserve solutions of (6). It is standard 
to denote the action of v on solutions u(x,t) of (6) by 

p(expev)t/(x,f) = a(x,t;e)u(ai(x,t;e),a2(x,t;e)) (8) 

for some functions a, a\, and «2, where e is the parameter of the group, a is referred 
to as the multiplier, and a\ and a 2 are changes of variables of the symmetry. For the 
applications we have in mind, e and a are of crucial importance, e will play the role 
of the transform parameter of the Fourier or Laplace transform and a will usually 
be the Fourier or Laplace transform of the transition density. Following [19], we 
assume that (6) has a fundamental solution p(t ,x,y). For this paper, it suffices to 
recall that we can express a solution u(x,t) of the PDE (6) subject to the initial 
condition u(x,0) = f(x) in the form 



u {x,t) = / f(y)p(t,x,y)dy, 
Jn 



(9) 



where p(t,x,y) is a fundamental solution of (6). The key idea of the transform 
method is to connect (8) and (9). Now consider a stationary, i.e. a time-independent 
solution, say uq(x). Of course, (8) yields 

p (expev) uq(x) = o(x,t;e)uo (a\(x,t;e)) , 

which also solves the initial value problem. We now set t = and use (8) and (9) to 
obtain 

/ a(y,Q,e)uo(ai(y,Q,;e)) p(t,x,y)dy = a (x,t;e)u(j(ai(x,t;e)) . (10) 
Jn 

Since a, uq, and a\ are known functions, we have a family of integral equations for 
p(t,x,y). To illustrate this idea using an example, we consider the one-dimensional 
heat equation 

1 9 

We will show that if u(x,t) solves (11), then for e sufficiently small, so does 

et 2 ze 



u(t,z) =exp<j -pr \ u{z-t£,t) 



Taking uq = 1, (10) gives 
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{y£ 1 f t£ 2 x£ 

-p- j p(t,x,y)dy = exp |^-2 - -j 

Setting a = —\, we get 

f" { a 2 g 2 t ] 

/ exp{ay}p(t,x,y)dy = exp i +ax\ . (12) 

We recognize that (12) is the moment generating function of the Gaussian distri- 
bution, so p{t,x,y) is the Gaussian density with mean x and variance g 2 t. We alert 
the reader to the fact that e plays the role of the transform parameter and a corre- 
sponds to the moment generating function. Finally, we recall a remark from [17], 
namely the fact that Laplace and Fourier transforms can be readily obtained through 
Lie algebra computations, which suggests a deep relationship between Lie sym- 
metry analysis and harmonic analysis. Lastly, we remark that in order to apply the 
approach, we require the PDE (6) to have nontrivial symmetries. The approach de- 
veloped by Craddock and collaborators, see [17], [18], [19], [20], and [21], provides 
us with the following: A statement confirming if nontrivial symmetries exist and an 
expression stemming from (10), which one only needs to invert to obtain p(t ,x, y). 
We first present theoretical results, and then apply these to the case of the MMM. 
Now we discuss the question whether the PDE (6) has nontrivial symmetries, see 
[20], Proposition 2.1. 

Theorem \.Ify^2, then the PDE 

u, = bx^u xx + f(x)u x — g(x)u , x>Q,b>0 (13) 

has a nontrivial Lie symmetry group if and only if h satisfies one of the following 
families of drift equations 

bxh' -bh + hi 2 + 2bx 2 - r g(x) =2bAx 2 -y + B, (14) 

bxti-bh+ l -h 2 + 2bx 2 -y g {x) + (15) 



bxti - bh + \h 2 + 2bx 2 -~ i g{x 



with K=l (y-4)b 2 andh(x) =x 1 ~ r f(x). 



(16) 



For the case y = 2, a similar result was obtained in [20], Proposition 2. 1 . Regarding 
the first Ricatti equation, (14), the following result was described in [20], Theorem 
3.1: 

Theorem 2. Suppose 7^2 and h{x) =x i ~Yf(x) is a solution of the Ricatti equation 
bxh' -bh + hi 2 + 2bx 2 -^g{x) = 2bAx 2 - y + B . 
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Then the PDE (13) has a symmetry of the form 



e(,) "^^ P l ^(1+4*) J (17) 

ex p|^(^( ^T-l-^Wll (18) 



{i+4t-t)' I= y 



37' l+4« 



(19) 



where F'(x) = f(x)/x y and u is a solution of the respective PDE. That is, for £ 
sufficiently small, U £ is a solution of (13) whenever u is. Ifu(x,t) — uq (x) with uq an 
analytic, stationary solution there is a fundamental solution p(t,x,y) of (13) such 
that 

/ e,xv{-Xy 2 ~ y }uo{y)p{t,x,y)dy = U x (x,t). 
Jo 

Here Ui(x,t) = V i^-yf-x- Further, ifuQ = 1, then JJ° p(t,x,y)dy = 1. 

For the remaining two Ricatti equations, (15) and (16), we refer the reader to Theo- 
rems 2.5 and 2.8 in [17]. 

We would now like to illustrate how the method can be used. Consider a squared 
Bessel process of dimension S, where 5 > 2, 



dX, = Sdt + 2^ t dW, , 

where Xo = x > 0. The drift f(x) = 8 satisfies equation (14) with A = 0. Conse- 
quently, using Theorem 2 with A = and u(x,t) = 1, we obtain 

{Aex 1 _s 
-WT^i)Y l+Aety ^ 

where b = 2. Setting e = we obtain the Laplace transform 

Ux(x,t)=f vup{-Xy}p(t,x,y)dy 
Jo 

which is easily inverted to yield 

where v = § — 1 denotes the index of the squared Bessel process. Equation (20) 
shows the transition density of a squared Bessel process started at time in x for 
being at time t in y. Recall that I v denotes the modified Bessel function of the first 
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kind. This result, together with the real world pricing formula, (5), allows us to price 
a wide range of European style and path-dependent derivatives with payoffs of the 
type H = f(S* { ,S[ 2 ,. . . ,S*,), where d > 1 and ti,t2,---,t<i are given deterministic 
times. 

By exploiting the tractability of the underlying processes, Lie symmetry methods 
allow us to design efficient Monte Carlo algorithms, as the following example from 
[2] and [3] shows. We now consider the problem of pricing derivatives on realized 
variance. Here we define realized variance to be the quadratic variation of the log- 
index, and we formally compute the quadratic variation of the log-index in the form, 

Recall from Section 2 that Y = {Y, , t > 0} is a square-root process whose dynamics 
are given in equation (2). In particular, we focus on put options on volatility, where 
volatility is defined to be the square-root of realized variance. We remark that call 
options on volatility can be obtained via the put-call parity relation in Lemma 4. 1 in 
[2]. The real-world pricing formula (5) yields the following price for put options on 
volatility 

( K ~ yj T So T s , 



sf'E 



S j- 



4 ■ (21) 



For computing the expectation in (21) via Monte Carlo methods, one first needs 
to have access to the joint density of (Sj*,Jq and subsequently perform the 
Monte Carlo simulation. Before presenting the relevant result, we recall that Sj* = 
SjOCpYT, i.e. it suffices to have access to the joint distribution of (Yr,jJ *). We 
remark that if we have access to the Laplace transform of (Yj, f i.e. 

E(exp(-AF T - M /f)), (22) 

then we have, in principle, solved the problem. From the point of view of imple- 
mentation though, inverting a two-dimensional Laplace transform numerically is 
expensive. The following result from [20], see Corollaries 5.8 - 5.9, goes further: In 
fact the fundamental solution corresponds to inverting the expression in (22) with 
respect to A, which significantly reduces the computational complexity. 

Lemma 1. The joint Laplace transform ofYj and Jq j- is given by 
_r(3/2 + v/2)„ , 



fix exp [7] \T 



r(v + l) r \ \ tanh(77772) 
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1 „_, a 2,„,„M„ (P 



where a = T] (l+coth(^)) +A, j3 = jtfrgy v = 2 y \ + 2 M> andM, %r (z) denotes 

the Whittaker function of the first kind. In [20], the inverse with respect to X was 
already performed explicitly and is given as 

p{T,x,y) - 



sinh (77T/2) \x; 

A"( T+ '- y -^um)H^m)- <23) 

Consequently, to recover the joint density of (Yt,Jq *), one only needs to invert 
a one-dimensional Laplace transform. For further details, we refer the interested 
reader to [3]. By gaining access to the relevant joint densities, this example demon- 
strates that Lie symmetry methods allow us to design efficient Monte Carlo algo- 
rithms for challenging finance problems. 



4 Wishart Processes 

Very tractable and highly relevant to finance are models that generalize the pre- 
viously mentioned MMM. Along these lines, in this section we discuss Wishart 
processes with a view towards exact simulation. As demonstrated in [13], Wishart 
processes turn out to be the multidimensional extensions of squared Bessel pro- 
cesses. However, they also turn out to be affine, see [33], and [34]. Prior to the latter 
two contributions, the literature was focused on affine processes taking values in 
the Euclidean space, see e.g. [27], and [28]. Subsequently, matrix-valued affine pro- 
cesses were studied, see e.g. [22], and [35]. Since [33], and [34], it has been more 
widely known that Wishart processes are analytically tractable, since their charac- 
teristic function is available in closed form; see also [31]. In this section, we exploit 
this fact when we discuss exact simulation of Wishart processes. 

Firstly, we fix notation and present an existence result. Wishart processes are 
S 1 ! or S~Jj valued, i.e. they assume values in the set of positive definite or positive 
semidefinite matrices, respectively. This makes them natural candidates for the mod- 
eling of covariance matrices, as noted in [33]. Starting with [33] and [34], there is 
now a substantial body of literature applying Wishart processes to problems in fi- 
nance, see [14], [15], [23], [24], [25], [26], and [32]. In the current paper we study 
Wishart processes in a pure diffusion setting. For completeness, we mention that 
matrix valued processes incorporating jumps have been studied, see e.g. in [5], and 
[46]. These processes are all contained in the affine framework introduced in [22], 
where we direct the reader interested in affine matrix valued processes. 
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In the following, we introduce the Wishart process as described in the work of 
Grasselli and collaborators; see [25] and [35]. For x G St, we introduce the 
valued Wishart process X x = X = {X t ,t > 0}, which satisfies the SDE 

dX, = (aa T a + bX t +X t b T ^dt+(yT t dW t a + a T dWjVx^j , (24) 

where a > 0, b G ^t&, a G ^Md- Here denotes the set of d x d matrices taking 
values in 3i. An obvious question to ask is whether equation (24) admits a solution, 
and, furthermore, if such a solution is unique and strong. For results on weak solu- 
tions we refer the reader to [22] , and for results on strong so lutions to [51]. We now 
present a summary of results, which in this form also appeared in [1]; see Theorem 
1 in [1]. 

Theorem 3. Assume thatx G St, and (X >d—\, then equation (24) admits a unique 
weak solution. If x G St and a > d + 1, then this solution is strong. 

In this paper, we are interested in exact simulation schemes to be used in Monte 
Carlo methods. Hence weak solutions suffice for our purposes and we assume that 
a > d— 1, so that the weak solution is unique. As in [1], we use WISd(x, a,b,a) to 
denote a Wishart process and WISd{x, a,b,a;t) for the value of the process at the 
time point t. 

We begin with the study of some special cases, which includes an extension of 
the MMM to the multidimensional case. We use B t to denote an « x d Brownian 
motion and set 

X, = BjB, . (25) 
Then it can be shown that X = {X t , t > 0} satisfies the SDE 

dX, = nl d dt + VxldW, + dWj VT, , 

where W t is a d x d Brownian motion, and I e i denotes the d x d identity matrix. This 
corresponds to the case where we set 

a = I c i , b = , a = n . 

We now provide the analogous scalar result, showing that Wishart processes gener- 
alize squared Bessel processes: Let 8 G J/ , and set 

k=\ 

Now we set 

5 

X t =Y,(W, k + w k ) 2 . (26) 

k=\ 

Then X can be shown to satisfy the SDE 



dX, = 8dt + 2VX,dB t7 
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where B = {B, , t > 0} is a scalar Brownian motion. This shows that (25) is the 
generalization of (26). Furthermore, it is also clear how to simulate (25). 

Next, we illustrate how Wishart processes can be used to extend the MMM from 
Section 2. We recall some results pertaining to matrix-valued random variables, see 
e.g. [36], and [52]. We introduce some auxialiary notation. We denote by ^# m ,„(9J) 
the set of all m x n matrices with entries in SR. Next, we present a one-to-one rela- 
tionship between vectors and matrices. 

Definition 1. Let A G ^# mj „ (3i) with columns a, € 9t m , i = 1 , . . . , n, and define the 
function vec : Jt m .n^) -> SR m " via 




We can now define the matrix variate normal distribution. 

Definition 2. A p x n random matrix is said to have a matrix variate normal distri- 
bution with mean M g ^# p ,„(5R) and covariance E®W, where E e S"p, W £ ^„ + , 
if vec(X J ) ~ ^4 / pn (vec(M J ),E <g> f ), where jV pn denotes the multivariate normal 
distribution on 5R P " with mean vec(M T ) and covariance E £§) f. We will use the 
notation X ~ jY pfl (M, E&W). 

Next, we introduce the Wishart distribution, which we link in the subsequent theo- 
rem to the normal distribution. 

Definition 3. A p x /^-random matrix X in is said to have a noncentral Wishart 
distribution with parameters p 6 jY , n> p, E S S^t and & e if its proba- 

bility density function is of the form 

fx(s) 

2^r p fydet(E)?y l etr^(0 + E- l S) 

det{S)^"- p - l '> Fi (^- l -0E- l S 

where S £ ,5^ p and qF\ is the matrix-valued hypergeometric function, see [36], and 
[52] for a definition. We write 

X ~ W p (n,E,&). 

Before stating the next result, recall that scalar non-central chi-squared random vari- 
ables of integer degrees of freedom, can be constructed via sums of normal random 
variables; see e.g. [41]. The following result presents the matrix variate analogy. 

Theorem 4. LetX ~ Jf P) „(M,E®I n ), n G {p,p+l,...}. Then 
XX T ~W p (n,E,E- 1 MM T ). 
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5 Bivariate MMM 



Theorem 4 is now employed to extend the MMM to a bivariate case. We consider 
exchange rate options, and follow the ideas from [37]. The GOP denominated in 
units of the domestic currency is denoted by S a , and the GOP denominated in the 

foreign currency by S . An exchange rate at time t can be expressed in terms of a 

^ s" 
ratio of two GOP denominations. Then one would pay at time t , ■+ units of currency 

a to obtain one unit of the foreign currency b. As the domestic currency is indexed 
by a, the price of, say, a call option with maturity T on the exchange rate can be 
expressed via the real world pricing formula (5) as: 



S a E 



\ 



(27) 



We now discuss a bivariate extension of the MMM from Section 2, which is still 
tractable, as we can employ the non-central Wishart distribution to compute (27). 
For k £ {a,b}, we set 

ak _ aO,kpk 

O, — O, Of , 

where sf' = expjr^f}, S^' = 1, so S°' k denotes the savings account in currency 
k, which for simplicity is assumed to be a deterministic exponential function of 
time. As for the stylized MMM, we model the discounted GOP, S*, denominated in 
units of the A:th savings account, S®' k , as a time-changed squared Bessel process of 
dimension four. We introduce the 2x4 matrix process X = {X t , t > 0} via 




The processes W^'', i= 1, ... ,4, denote independent Brownian motions, subject to 
the deterministic time-change 

<p 1 (0-^ r (exp{ J7 1 f }-l) = i| o 'a>,, 

c.f. Section 2. Similarly, also W l f , i= 1, ... ,4, denote independent Brownian mo- 
tions, subject to the deterministic time change 

9 2 W = ^(ex P {r ? 2 /}-l) = 5/ '«>- 

Now, consider the process Y = {Y t , t > 0}, which assumes values in St, and is 
given by 
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Y, :=X,Xj ,t >0, 



which yields 



Y t = 



- w 



Li 



ULiLU {Km 



■w 



l,J 



- 11" 



We set 



and 



i,2 



sf = y r 2 - 2 , 



so we use the diagonal elements of 7, to model the GOP in different currency de- 
nominations. Next, we introduce the following dependence structure: The Brownian 
motions W 1 and W' 2 , i= 1 , ... ,4, covary as follows, 



(28) 



where —1 < p < 1. The specification (28) allows us to employ the non-central 
Wishart distribution; we work through this example in detail, as it illustrates how to 
extend the stylized MMM to allow for a non-trivial dependence structure, but still 
exploit the tractability of the Wishart distribution. We recall that vec(Xj) stacks the 
two columns of Xj, hence 



vec{X~l) 



VV mUl 



,1,1 



4.1 

L,2' 



It is easily seen that the mean matrix M of vec(Xj) satisfies 

,1,1 ' 



,1.2 



w 



4.2 



(29) 
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(30) 



and the covariance matrix of vec{Xj. ) is given by 

r . _ rz u / 4 E i < 2 i 4 

where E is a 2 x 2 matrix with £ u = cp 1 (r), E 2 - 2 = (p 2 (T), and 

We remark that assuming — 1 < p < 1 results in E being positive definite. It now 
immediately follows from Theorem 4 that 



X T X\- ^W 2 [^,E,E~ l MM' 
where M and E are given in equations (29) and (30), respectively. Recall that we set 



Y, 


= x t xj 


sf 


_yl.l 


s» 


y2,2 



hence we can compute (27) using 

E(f(Y T )), 

where / : — > 5? is given by 



/(y) 



/ exp{/-ir}v 1 -' ^ 
\exp{r27"}v 22 



exp{riJ}y 



l.i 



for y S S2, and y''', / = 1,2, are the diagonal elements of y, and the probability 
density function of Yj is given in Definition 3. 

We now discuss further exact simulation schemes for Wishart processes, where 
we rely on [1] and [6]. For integer valued parameters a in (24), we have the fol- 
lowing exact simulation scheme, which generalizes a well-known result from the 
scalar case, linking Ornstein-Uhlenbeck and square-root processes. In particular, 
this lemma shows that, in principle, certain square-root processes can be simulated 
using Ornstein-Uhlenbeck processes. 

Lemma 2. Let A > 0, Q > 0, and define the SDEs 

dX; = -AX\dt + QdWi , 

for i = 1 , . . . , p\ where £ jY , W l , W 2 , . . .,W& are independent Brownian motions. 
Then 
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i=l 

is a square-root process of dimension j3, whose dynamics are characterised by an 
SDE 

dZ, = (PQ 2 - 2AZ t )dt + 2Qy/ZdB, , 
where B is a resulting Brownian motion. 

Proof. The proof follows immediately from the Ito-formula. □ 

This result is easily extended to the Wishart case, for integer valued a, see Section 
1.2.2 in [6]. We define 

P 

V t = ZVlt> (31) 

k=l 

where 

dX kJ = AX kJ dt + Q T dW kJ ,k=l,...,P, (32) 

where A G J£ d , X, G 3i d , Q e Ji d , W k G 9i d , so that V t G Jtd- The following lemma 
gives the dynamics of V = {V t , t > 0}. 

Lemma 3. Assume that V t is given by equation ( 31 ), where X, satisfies equation 
(32). Then 

dV t = ^Q T Q+AV, + V t A T ^dt + ^/V~ t dW t Q + Q T dWjy / V), 

where W = {W t ,t>0}isadxd matrix valued Brownian motion that is determined 
by 

P 

W~ t dW t =Y,Xk,tdWl k . 

k=l 

Finally, we remind the reader that vector-valued Ornstein-Uhlenbeck processes can 
be simulated exactly, see e.g. Chapter 2 in [55]. 

For the general case, we refer the reader to [1]. In that paper, a remarkable split- 
ting property of the infinitesimal generator of the Wishart process was employed to 
come up with an exact simulation scheme for Wishart processes without any restric- 
tion on the parameters. Furthermore, in [1] higher-order discretization schemes for 
Wishart processes and second-order schemes for general affine diffusions on posi- 
tive semidefinite matrices were presented. These results emphasize that Wishart pro- 
cesses are suitable candidates for financial models, since exact simulation schemes 
are readily available. 



6 Conclusion 



In this paper, we discussed classes of stochastic processes for which exact simu- 
lation schemes are available. In the one-dimensional case, our first theorem gives 
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access to explicit transition densities via Lie symmetry group results. In the mul- 
tidimensional case the probability law of Wishart processes is described explicitly. 
When considering applications in finance, one needs a framework that can accom- 
modate these processes as asset prices, in particular, when they generate strict local 
martingales. We demonstrated that the benchmark approach is a suitable framework 
for these processes and allows to systematically exploit the tractability of the mod- 
els described. For long dated contracts in finance, insurance and for pensions the 
accuracy of the proposed simulation methods is extremely important. 
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